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ABSTRACT 

Context. The afterglows and host galaxies of long gamma-ray bursts (GRBs) offer unique opportunities to study star-forming galaxies in the high-z 
Universe. Until recently, however, the information inferred from GRB follow-up observations was mostly limited to optically bright afterglows, 
biasing all demographic studies against sight-lines that contain large amounts of dust. 

Aims. Here we present afterglow and host observations for a sample of bursts that are exemplary of previously missed ones because of high visual 
extinction (A™^ > 1 mag) along the sight-line. This facilitates an investigation of the properties, geometry and location of the absorbing dust of 
these poorly-explored host galaxies, and a comparison to hosts from optically-selected samples. 

Methods. This work is based on GROND optical/NlR and Swift/XKT X-ray observations of the afterglows, and multi-color imaging for eight 
GRB hosts. The afterglow and galaxy spectral energy distributions yield detailed insight into physical properties such as the dust and metal content 
along the GRB sight-line as well as galaxy-integrated characteristics like the host's stellar mass, luminosity, color-excess and star-formation rate. 
Results. For the eight afterglows considered in this study we report for the first time the redshift of GRB 081 109 (z = 0.9787 ± 0.0005), and the 
visual extinction towards GRBs 081 109 (A°'^'' = 3.4!°* mag) and 100621A (A°^^ = 3.8 ± 0.2 mag), which are among the largest ever derived for 
GRB afterglows. Combined with non-extinguished GRBs, there is a strong anti-correlation between the afterglow's metals-to-dust ratio and visual 
extinction. 

The hosts of the dustiest afterglows are diverse in their properties, but on average redder ({(R - A')ab> ~ 1-6 mag), more luminous ({L) ~ 0.9 L*) 
and massive ((log M, [Mq]) ~ 9.8) than the hosts of optically-bright events. We hence probe a different galaxy population, suggesting that previous 
host samples miss most of the massive, chemically-evolved and metal-rich members. This also indicates that the dust along the sight-line is often 
related to host properties, and thus probably located in the diffuse ISM or interstellar clouds and not in the immediate GRB environment. Some 
of the hosts in our sample, are blue, young or of small stellar mass illustrating that even apparently non-extinguished galaxies possess very dusty 
sight-lines due to a patchy dust distribution. 

Conclusions. The afterglows and host galaxies of the dustiest GRBs provide evidence for a complex dust geometry in star-forming galaxies. In 
addition, they establish a population of luminous, massive and correspondingly chemically-evolved GRB hosts. This suggests that GRBs trace the 
global star-formation rate better than studies based on optically-selected host samples indicate, and the previously-claimed deficiency of high-mass 
host galaxies was at least partially a selection effect. 
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1. Introduction 



Long gamm a-ray bursts (GRBs, see e.g., [Zhang |[2007l |Gehrels| 
|et al.|2009 for reviews) are linked to core-collapse supernovae 
and hence star-formation via the death of massive stars (e.g., 
[Galama et aL][T998l [Hjorth et ai;i[2003] l. At high redshifts, 
where a significant fraction of the star-formation is thought to 
be dust-obscured (e.g., [Adelberger & Ste idel 200()| [Chapman 



|et aH2005[ ), GRBs and their host galaxies offer an independent 
track towards a better understanding and full census of the star- 
formation in the early Universe: GRBs, having luminous emis- 
sion in a simple power-law spectrum provide the ideal back- 
ground light to illuminate dust-enshrouded star-forming regions 
which would otherwise remain unexplored, while at the same 
time pinpointing their host galaxies. 

However, the extent to which GRB hosts provide an unbiased 
picture of the formation of high-mass stars, and whether they 

* Based on observations made with GROND at the MPl/ESO 2.2 m 
telescope and with telescopes at the European Southern Observatory at 
LaSilla/Paranal, Chile under program 086.A-0533 and obtained from 
the ESO/ST-ECF Science Archive Facility from programs 177.A-0591 
and078.D-0416. 
** Visiting Astronomer at MPE 



preferentially occur in low-metallicity environments remains a 
much debated issue (e.g., Le Floc'h et al. 2003; Fynbo et al. 
|2003l ITanvir et al.|[2504l |Fruchter et aTpOOe, Kocevski~eraL 



2009 1. In single progenitor models, metal-poor stars are favored 



by theory ( Woosley 


1993 


MacFadyen & Woosley 1999 1, as they 


would in principle 


5e ab 


e to keep more angular momentum at 



the time of stellar collapse due to smaller wind pressures and 
losses throughout their evolution (e.g., lYoon & Langer 2005 



|Mokiem et al.|2007) l. However, binary progenitor channels could 
also play an important role in the formation of long GRBs (e.g.. 



Fryer et al. 1999 ), having somewhat relaxed metallicity con- 
straints relative to single star progenitors ( [Fryer et al.||2007) . 
Observations of GRB hosts are hence not only important in a 
cosmological context, but provide relevant clues to the exact na- 
ture of GRB progenitors. 

A fundamental limit of hitherto available GRB host galaxy 
samples is the incompleteness which arises from the non- 
detection of the optical afterglow of a GRB (e.g., Groot et al. 
|1998[ Fynbo et al.1'2001 ). These optically dark bursts could be 
caused by either high-redshift (e.g., Grei ner et al.|2009| Tanvir 
etar2009; S alvaterra et"ar.| 2009, Cucchi ^aTt al'|2011|l large 



columns of dust (e.g., [Klose et al.|2000t|20()3flT^nvir et al.|2008[ 



2 
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|Perley et al.|2011bl l or an intrinsically fainter optical afterglow 
as compared to the extrapolation of X-ray data when using syn- 
chrotron emission theory, i.e., a decoupled optical/X-ray after- 
glow light-curve (e.g., |Panaitescu et al. 2006; Ghiselli ni et al.| 
2009 |Nardini et al. 2010V New afterglow samples became avail- 



able through the recent advent of dedicated afterglow follow-up 
campaigns on medium-to-large aperture telescopes (e.g., |Fynbo| 
[eral',2009a, ,Cenko et al. 2009; Greineret al.|2011| l . Th ese new 
afterglow samples reach completeness levels of ~90% ( Greiner 
et al.|20lT I and settled the dark-burst issue: Around three quar 



ters of dark bursts are the result of a dusty afterglow line of sight 
(e.g., Perley et al.|2 009 ; Gre iner et al.|201 1 ). Accurate positions 
from afterglow observations are necessary to unambiguously as- 
sociate galaxies to GRBs. The lack of optical/NIR afterglows for 
dark GRBs therefore implies an inherent bias against the associ- 
ated host galaxies. 

The available host population is not as complete as the most 
recent afterglow samples. Instead, it is largely selected from op- 
tically bright afterglows and shows a prevalence of young and 
vigorously star-forming galaxies wit h sub-L* luminosities and 
masses around 10^ Mq (e.g.. Blo om et al.|l998f,Le Floc'h et al 



2003 



Christensen et al. 2004 ; Fruchter et al. 2006 Savaglio et al. 



20091 referenced as SGL09, hereafter). However, it is an open 
question whether this is a physical consequence of GRBs prefer- 
ring low-metallicity environments, or merely a selection effect: 
Host galaxies of dark GRBs were typically not identified, and 
hence are under-represented in the available host sample. 

Whether the physical characteristics of hosts of optically 
dark and bright GRBs are distinct is also the subject of discus- 



sion. Previous sample studies (e.g., Berger et al.|2003 Le Floc'h 



et al.|2003 Perley et al. 20091 have not revealed strong evidence 



of a significant difiference. A handful of single dark GRBs were 



ities and stellar masses over 10" Mq (e.g.. 


Levan et al.||2006 


Berger et al.|2007||Levesque et al.|2010 


IHashimoto et al.|2010 


Ktipcii Yolda§ et al. 2010 Chen et al. 


20101. Recently Perley 


et al. (2010a |2011a|l suggested, that the 


general galaxy popula- 



tion hosting dark bursts is redder and more luminous and hence 
suggestive of a higher metallicity than those selected via opti- 
cally bright afterglows. 

In this paper, we study the nature of GRB hosts that previ- 
ously escaped detection due to the dust bias, and are hence ex- 
emplary of those missing from demographic studies. We avail of 
dedicated GRB afterglow campaigns with high completeness to 
preselect the GRB hosts for this study. These afterglow data do 
not only provide accurate positions for host identifications, but 
for the first time allow us to directly select dust extinguished 
(and not only optically faint) GRBs via well-sampled broad- 
band (NIR to X-ray) afterglow observations. 

After selecting afterglows with visual extinctions Ay™ ex- 
ceeding unity, we search for the associated hosts with the 
Gamma-Ray Optical and Near-Infrared Detector (GROND, 
Greiner et al., 200 8) as well as the ESO Very Large and New 
Technology Telescopes (VLT and NTT, respectively) and the 
Ultra- Violet Optical Tel escope (UVOT, [Roming et al.|2005] l on- 



board the Swift satellite ( [Gehrels et al.|2004| l. 

The obtained data allow us to study the physical properties 
of the hosts of high-A^"^^ GRBs in detail, and to investigate the 
bias against dust in GRB host samples. As an ultimate conse- 
quence they address the role of GRB host galaxies as tracers 
of galaxy formation and evolution. Furthermore, they link after- 
glow diagnostics, i.e., detailed information about a single sight 
line, to host-integrated properties, and in this combination di- 



10" 



I I Previous Afterglows 
This Sample 



I: 



Afterglow Ay [mag] 

Fig. 1. Histogram of visual extinction values from previous af- 
terglows and those selected for this study. We note that five of 
the afterglows in this work (GRBs 070802, 080605, 080607, 
080805, 090926B) were already discussed in previous sample 
studies, but are only included in one of the histograms for clar- 
ity. 



rectly probe the nature of dust and its properties in high-redshift, 
star-forming galaxies. 

Throughout this work, we adopt the convention that the flux 
density of the afterglow Fy(v, t) can be described as Fv(y, i) oc 
y-l^r", and concordance (Qm = 0.27, Qa = 0.73, Hq = 
71 km/s/Mpc) ACDM cosmology. All errors are given at Icr con- 
fidence unless indicated otherwise. All magnitudes and colors 
are given or converted into the AB system. 



2. Sample Selection 

The host galaxy sample presented in this work is based on a 
direct measurement of large visual extinction along the GRB 
line of sight {A^™ > 1 mag) from multi-color (NIR to X- 
ray) afterglow observations. Specifically, eight GRB afterglows 
(GRBs 070306, 070802, 080605, 080607, 080805, 081109, 
090926B and 100621A) fulfill the selection criterion and define 
our initial host sample. Our host sample is a direct result of af- 
terglow observations. The selection itself is hence not limited by 
galaxy brightness, nor introduces a bias towards luminous galax- 
ies. Afterglow measurements for the initial selection have been 
obtained from the literature or by analyzing photometric opti- 
cal/NIR data from the GROND archive. In the latter case, they 
are detailed in Sect. 13.11 

Our eight GRBs have a median redshift of {zav) — 1-5- 
This is significantly lower than the published mean of Swift 
GRBs with measured redshifts {{zs 



,if,) - 1.9, Fynbo et al. 
|2009 f|, but larger than the one of the previous host sample 
((zsgloq) = 0.96 from SGL09), which includes a large number 
of pre-5w//r events. Within the selected eight afterglows, four of 
them show a clear 2175A dust feature (GRBs 070802, 080605, 
080607, 080805), which is with flie exce ption of GRB 080603 A 
( [Guidorzi et al.|[20lT Kann et aL]|201 1)) the full sample where 
a sig nificant detection o f this feature has been reported to date 
(e.g., Zafar et al. 2011 1. In fact, the 2175 A feature could not 



' Including updates from http://www.mpe.mp g.de/~l jcg/grbgen.html 
and http://www.raunvis.hi.is/~pja/GRBsample.html 
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have been detected in the rest of the sample. This arises from 
the combination of a large dust column, and not deep and rapid 
enough follow-up in the case of GRBs 070306 and 090926B, or 
no observational wavelength coverage at 2175 A x (1 +z) (GRBs 
081109, 100621 A). 

While observationally challenging, the requirement of a dust 
measurement has several obvious advantages over a selection 



based on the optical-to-X-ray flux ratio (j6ox, see e.g., Jakobsson 



let al.|2004| |Rol et al.|2005) |van der Horst et al.|2009| l. Most im- 
portantly, the selection is the result of a measurement rather than 
an extrapolation, and provides a clean selection of dusty GRBs: 
Our afterglows are chosen accord to their visual extinction, in- 
stead of only their optical faintness. The latter events could of 
course also be at high-z or have different emission mechanisms 
in the optical and X-ray regime. 

We note, that our selection is still somewhat model- 
dependent, specifically that the afterglow emits synchrotron ra- 
diation in the optical/NIR and X-ray regime. Despite the lack 
of conclusive alternatives to the standard synchrotron-fireball 
model, there are still puzzling features in well-sampled multi- 
color light curves, such as chromatic]^ breaks, optical and/or X 



ray flares or plateaus (e.g., Panaitescu et al. 
[2009'; Kruhl er et al.|[2509l [Gates et al.||201 1 



2006 Evans et al. 



^ ^ ^ ^ also illustrated in 

Fig.|3]l. The apparently decoupled optical and X-ray light curve 
for some bursts results in a strong dependence of /3ox on the time 
of the observation (which in fact is directly observed in some 
cases, see e.g., Filgas et al.|2011[ ), and hence the dark burst def- 
inition depends also on observational constraints and not only 
physical properties. 

In addition, the required optical and/or NIR afterglow detec- 
tion yields the GRB position to sub-arcsec accuracy, and hence 
negligible chance-coincidence probabilities for field galaxies, 
which is particularly important for small sample sizes as in this 
work. For all events a spectroscopic redshift is available from 
the literature, or could be obtained via host galaxy spectroscopy. 
This enables quantitative studies, such as the comparison be- 
tween specific sight-lines against host-integrated properties, as 
well as an investigation of the relation between the dusty hosts 
to the host sample of SGL09. 

Due to the inherent difficulties to accurately localize high- 
afterglows, most previous afterglow samples are biased 
towards small visual extinctions. This is illustrated in Fig. [T] 
which shows the visual extinction of GRBs selected for this work 
as compared to previously compiled Ay™ values. The latter 
were taken from,Kann et al.|(2006,,20Tof;ISchady et al.[ ( [20T0) i; 



/I GRB 

Ay 



[Greiner et al.1 ( [2()l l| i and [Zafar et al.[ ( [20"TT 



3. Observations 

3.1. Afterglows 

Optical and near-infrared measurements of the afterglows of 
GRBs 070306, 070802, 080605, 080607, 080805 and 090926B 
or results thereof are taken from 'Jaunsen et al.| ([2008[l, Kruhler[ 



etaL] ([200 8 ), Perl ey et al.]p)T lb), Grei ner et al.| ( [2011[ ), and 
Zafar et al.| (|201 1 ), respectively. GROND observations of the af- 
terglows of GRB 081109 and GRB 100621 A are not presented 
elsewhere and are shortly described in the following. 



^ Chromatic light-curve breaks are breaks in the X-ray regime not 
associated with contemporanous breaks in the optical/NIR bands, and 
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Fig. 2. GROND //-band fight curve of the afterglow of 
GRB 081109. Highlighted with a grey shaded area and labeled 
with I is the time interval which has been used to extract a si- 
multaneous, host subtracted afterglow SED from the GROND 
and SwiftlXWY data. 



3.1.1. GRB 081109 



Swift triggered on GRB 081 109 C Immler et al.|2 008'), and X-ray 
and NIR measurements of the afterglow were rapidly reported 
by [Beardmore et al.[p008] l and [D'Avanzo et al [ ( l2U08] l. GROND 
observations in seven optical/NIR filters (g'r'i'z JHKg) simulta- 
neously, started 17.1 hr after the GRB trigger (Clemens et al. 
|2008b| l and a preliminary analysis of the spectral energy distri- 
bution (SED) already revealed significant reddening of the op- 



tical/NIR afterglow ( [Clemens et a l. 2008a). GROND continued 
to observe the transient at 2, 3, 6 and 378 days after the trig- 
ger, where the host brightness was derived from the last epoch. 
GROND afterglow and host measurements are given in Table [T] 
and Tables[3]and|4] respectively. A possible host galaxy was also 
reported in the UVOT white filter ( |Kuin & Immler|2008| ). High- 
energy prompt and afterglow data, early NIR imaging and a light 
curve and X-ray spectral analysis of this burst are presented by 
[Jin et al.[ ( |2009] l. 

3.1.2. GRB 100621A 



GROND reacted immediately (" Updike et"aL][20 10l) to the Swift 
trigger of GRB 100621 A (Ukw atta et al.[ [2010b) taking the 
first images at 230 s after the burst. Simultaneous imaging in 
g'r'i'z'JHKs continued for 3.05 h, and was resumed on night 2, 
4 and 10 after the burst. Analysis of the Swift/X-Tay data and 
further detections of the NIR afterglow were reported by |Stratta 
[etaL] ( [20T0| ) and [Naito etall (2010), respectively. Early GROND 
data already revealed evidence for substantial reddening and host 
emission dominating the flux in the bluest filters, which was 
also seen by Swift fVWOT ( [Gates & Ukwatta[[20T0| . Afterglow 
and host measurements are again shown in Tables [T] |3] and |4] 
The temporal evolution of the optical/NIR afterglow is complex 
with a very steep increase in brightness of around 1.5 mag in 
the J band from 3.5 to 4.5 ks after the trigger. The light curve 



vice versa. 



is hence very similar the one of GRB 081029 (Nardini et al. 
1201 1), where an analogous behavior could be associated with 
the intrinsic properties of the GRB and not to changes in the in- 
tervening dust content. For completeness, the /-band light curve 
of GRB 100621 A is shown in Fig.|3] but its detailed modeling 
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Fig. 3. GROND 7-band light curve of the afterglow of 
GRB 100621 A. Highlighted with grey shaded areas and labeled 
with I and 11 are the time intervals which have been used to ex- 
tract a simultaneous, host subtracted afterglow SED from the 
GROND and Swift/XRT data. 



and interpretation is beyond the scope of this paper, and will be 
discussed in a future work. 



3.2. Hosts 

Once the sample based on afterglow host extinction was defined, 
late follow-up observations were initiated first with GROND, 
and in case of non-detections in individual filters, were con- 
tinued with telescopes of successively increasing aperture size, 
specifically with EFOSC/SOFI at the NTT (4m class) and 
FORS2/HAWKI at the VLT (8m class). In one case without 
published redshift (GRB 081109), the photometric imaging was 
complemented by low-resolution spectroscopy with FORS2. 
Public VLT data for known hosts (GRBs 070306 and 070802) 
were obtained from the ESO archive. Ground-based data were 
complemented by Swift /WOT imaging for GRBs 081109 and 
100621A. In th e case of GRB 080607 all host measurements 
were taken from |Chen etaL| ( [20T0l[20TT] ). 



4. Data reduction and SED fitting 

4.1. Swift/XRT & UVOTdata 

X-ray data have been retrieved from the HEASARC archiv^ 
and reduced with the xrtpipeline. Spectra have been grouped to 
yield a minimum of 20 counts per bin, while the light curve was 
taken from the 5w!/f/XRT light curve repository dEvans et al. 
2007] [20091). Swiftl\]N (JY photometry has been obtained follow- 
ing Poole et al. (2008 1 and is provided in Table |2] 



4.2. Ground-based optical/NIR ptiotometry 

All optical (GROND, EFOSC, FORSl/2) and near-infrared 
(GROND, SOFI, HAWK-I, ISAAC) imaging wa s reduced in a 
standard manner using pyraf/IRAF ( [Tody 1993| l similar to the 



procedure outlined in Kriihler et al. (2008 1. For afterglow and 



host photometry, point-spread function (PSF) fitting and aper- 
ture photometry was used, respectively. The aperture diameter 



for individual hosts ranges between typically 1.0" to 2.5", cor- 
responding to values between 1 .5 to 4 times the PSF FWHM. It 
has been chosen sufliciently large to include the largest fraction 
of host flux, and given the typical extent of these galaxies (< 1", 
i.e., < 8.5 kpc at z ~ 1.5), the fraction of missed low-surface 
brightness emission is very likely not a major contribution to the 
presented measurements. 

PSF and aperture photometry were then flux calibrated 
against GROND observations of SDSS fields (Abazajian et al. 
2009) taken immediately before or after the GRB field for the 
optical g'r'i'z' filters. BVRI photometry was obtained by cre- 
ating a set of 20 - 30 secondary standards from the GROND 
photometry of field stars and the color terms from Luptoij^ U - 
band photometry was tied to observations of Landolt standard 
stars ( Landoltj 19921 ) taken during the same night at different air- 
masses, which allowed a reliable correction of the atmospheric 
extinction to be applied. 

NIR photometry was measured against 20 - 60 point sources 
from the 2MASS catalog ( [Skrutskie et al.|2006| ) in the 10' x 10' 
field of view of GROND. The zeropoint for the NIR imagers 
with smaller field of views (in particular SOFI and ISAAC) was 
then derived using 3-10 secondary standards common in the 
GROND and SOFI/ISAAC/HAWK-I frames. The HAWK-I Y- 
band imaging was calibrated against the z' and J measurements 
from field stars using the color term 

Y = 0.05 + 0.463 x(z-J) + J 

where all magnitudes are in the AB system. This color term 
has been obtained from synthetic photometry of stars with the 



http://heasarc.gsfc.nasa.gov/cgi-bin/W3Browse/swift.pl 



templates of Pickles ( 1998 ) and Chabrier et al. (2000), and yields 
an rms residual scatter for individual stars of 0.07 mag. 

This procedure resulted in typical absolute accuracies of 2- 
5% for the optical {Utoz') filters and 4-8% in the NIR (YJHKs), 
which have been added in quadrature to the error introduced 
by photon noise. All data used in the analysis have been cor- 
rected for the expected Galactic foreground extinction according 
to |Schlegel e"taL] ( [T998] ) with Ry = 3.08. For the selected fields, 
this correction is small, i.e., Eg-v S 0.15 mag in all cases. The 
uncertainty of order of 10% in the Galactic foreground correc- 
tion is hence not going to affect the overall results of this work. 
Ground-based photometric measurements of the afterglows and 
hosts are shown in Tables [T] [3] and |4] respectively. 

4.3. Long-slit spectroscopy 

In addition to the photometric observations, the host of 
GRB 081109 was also observed spectroscopically with the VLT 
equipped with FORS2. In total 2 x 1200 s spectra were obtained 
with the grisms 300y and 300/ and a long-slit width of 1.6". 
Acquisition images were taken in the V and / filter. The spec- 
troscopic data were obtained at airmasses of ~ 1.1 and see- 
ing of 0.9", which results in a line-spread function of approx- 
imately 1 .9 nm at 570.0 nm. The data were reduced using stan- 
dard procedures in pyraf/IRAF, with the wavelength solution ob- 
tained against an HeHgCd arclamp exposures with 25 lines leav- 
ing residuals of around 0.07 nm rms. Flux-calibration was per- 
formed against the spectro-photometric standard BPM1627'^ 
The wavelength- and flux-calibrated spectrum was corrected for 
Galactic foreground extinction and renormalized to the available 
photometry. 

http://www.sdss.org/dr7/algorithms/sdssUBVRITransform.htmI 
^ http://www.eso.org/sci/observing/tools/standards/spectra/bpml6274.html 
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Table 1. Afterglow photometric measurements 



GRB 


At [ks] 


g' 


r' 




z' 


J 


H 




081109 


65 


>24.53 


>24.75 


23.6(4) 


22.8(3) 


21.24(17) 


19.78(11) 


18.83(08) 


100621 A 


7.6 


21.56(17) 


19.92(08) 


18.67(06) 


17.71(04) 


16.02(04) 


15.02(05) 


14.32(06) 



Notes: At is the mean time of the observation after the GRB trigger All magnitudes are in the AB system and uncorrected for Galactic foreground 
reddening. Values in brackets correspond to photometric errors in units of valid digits. Upper limits are given at 3cr confidence. 

Table 2. Swift/UYOT UV {uvw2 to u) photometric measurements of GRB hosts 



Host 


uvwl 


uvm2 


uvwl 


u 


GRB 081109 


>24.1 


23.6(3) 


23.4(3) 


>22.9 


GRB 100621 A 


22.31(04) 


22.23(06) 


22.20(07) 


21.95(06) 



Notes: All measurements in the AB system and uncorrected for the Galactic foreground reddening. Values in brackets correspond to photometric 
errors in in units of valid digits. Upper limits are given at 3cr confidence. 

Table 3. Optical (U to z') photometric measurements of GRB hosts 



Host 
GRB 070306 
GRB 070802 
GRB 080605 
GRB 080805 
GRB 081109 
GRB 090926B 
GRB 100621 A 



23.23(14) 
23.71(13) 
21.95(10) 



S 

22.90(09) 
23.15(07) 



23.07(07) 
23.31(07) 
21.86(06) 



25.7(2) 
22.85(06) 



23.08(09) 

22.82(07) 

22.74(07) 
22.96(06) 
21.48(06) 



23.00(09) 
25.20(09) 

25.5(2) 



22.81(13) 
25.5(3) 
22.81(08) 

25.7(4) 
22.01(08) 
22.92(12) 
21.15(06) 



21.96(09) 



22.86(17) 

22.76(11) 

21.99(09) 
22.44(10) 
21.46(06) 



Notes: All magnitudes in the AB system and uncorrected for the Galactic foreground reddening. Values in brackets correspond to photometric 
errors in units of valid digits. Upper limits are given at 3cr confidence. Data for the host of GRB 080607 were taken from jChen et al.|p0T0|[20TT] l, 
and are not shown in the table. 



4.4. SED fitting 
AAA . Afterglows 

For the afterglow SED analysis. X-ray (0.3 - 10 keV) and opti- 
cal/NIR data were fit together under the assumption that the un- 
derlying continuum emission is well represented by synchrotron 
emission (e g.,|Sari et al.|1998||Galama & Wijers|200T] |Schady| 
let al.|2007| [Rossi et al.|20nal i: 

Rest-frame soft X-ray photons are absorbed by metals, pre- 
dominantly a-chain elements, while UV over optical to NIR 
wavelengths are decreasingly affected by dust absorption. Good 
coverage above 1 keV combined with NIR observations allows 
for an accurate determination of the continuum, and hence good 
constraints on the dust abundance (represented by the absorp- 
tion in the rest-frame, optical V-band, A^**^), the extinction law 
and total metal content (converted into a Hydrogen-equivalent 



column density Nh,x assuming solar abundances from Anders 



|& Grevesse„1989) along the line of sight. Specifically, single 
and broken power-law continua were used, where in the latter 
case the two power-law slopes Pi and ySa were fixed to yield 
jSi -I- 0.5 = y32 as expected for the cooling break of synchrotron 
emission in the slow cooling regime (e.g.. Sari et al. 1998 
IGranot & Sari|2002| ). 



The NIR to X-ray SEDs were fit in X-spec ( |Arnaud| 1996) 1 us- 
ing extinction laws for the Milky Way (MW), an d Sma ll (SMC) 
and Large (LMC) Magellanic Clouds from [Pel] (jl992|l. For the 
GROND photometry, measurements from the time frame indi- 
cated in Fig. |2] and [3] were used, where we chose interval II for 
GRB 100621A due to a better signal-to-noise ratio. Data taken 
at interval I yield consistent results for A^"^® and Nh,x- In the 
X-ray regime, where spectral changes in the late evolution of an 
afterglow are typically very moderate or in most cases even com- 
pletely absent, the full XRT data set with a constant hardness ra- 



tio was used to create a time-averaged spectrum. The latter was 
then rescaled to the flux value at the time of the GROND ob- 
servations derived from fitting the XRT light-curve with simple 
afterglow models. The early steep decay of the X-ray light curve 
and epochs of flaring activity were excluded from the combined 
spectrum as well as the light-curve fitting. 

4.4.2. Hosts 

UV/optical/NIR photometry of the hosts of the selected GRBs 
were analyzed in a standard way using stellar population synthe- 
sis (SPS) techniques to convert luminositi es into stellar masses 



M, (e.g., [Bell et al.|2003[|llbert et al.|2010i ) within LePhar^Jn 



detail, 3 x 10" galaxy templates based on models from |Bruzual 
|& Charlot[ ( [20031 ) with a universal IMF ( |Chabrier|2003] l and dif- 
ferent ages, star formation histories, extinction laws, reddening 
values and metallicities were fit to the data. In addition, emission 
lines are taken into account by converting the de-reddened UV 
luminosity into a star formation rate, and hence line strengths 
of Ly-g, Hg, H p, [OH] and [OIII] following [Kennicutt| 
and Ilbert et al. (2009). In particular for vigorously star-forming 



galaxies such as GRB hosts, this effect is a significant contribu- 



tion to even broad-band photometry (e.g., Watson et al.[ 2011 
and reaches values of up to ~ 0.2 - 0.3 mag ( Ilbert et al.|2009 1. 
For a direct comparison with results pubUshed in the literature 
(e.g., [Fontana et al . 2006 ; Marchesini et a l. 2009 ; ' Ilbert et al." 
[20T0l SGL09), die attenuation law from |Calzetti et al.[ ( |2000| l 
derived for starburst galaxies is used, unless different reddening 
laws provide a better fit to the host data at 90% confidence. We 
caution that access to the rest-frame NIR, which is the best tracer 
of the stellar mass of the galaxy, is somewhat limited for part of 
the sample. However, all hosts are detected in at least one filter 



http://www.cfht.hawaii.edu/~arnouts/LEPHARE 



6 T. Kriihler et al.: The SEDs and Host Galaxies of the dustiest GRB afterglows 

Table 4. NIR (Y to K) photometric measurements of GRB hosts 



Host 


Y 


-/grond 


J 


^GROND 


H 


^GROND 


K 


GRB 070306 




21.9(4) 


21.62(08) 


21.5(4) 


21.20(12) 


> 


21.1 


21.38(10) 


GRB 070802 






24.5(3) 










23.4(3) 


GRB 080605 




21.9(2) 




22.3(3) 




> 


21.1 




GRB 080805 






23.6(2) 










23.1(2) 


GRB 081109 


21.63(08) 


21.40(17) 


21.37(06) 


21.5(4) 




> 


20.6 


21.05(08) 


GRB 090926B 




22.3(4) 


21.88(13) 


> 21.6 


21.9(3) 


> 


20.9 


21.44(19) 


GRB 100621 A 


21.10(06) 


21.22(10) 


21.43(06) 


21.18(14) 




> 


21.1 


21.23(11) 



Notes: All measurements are given in the AB system and are uncorrected for the Galactic foreground reddening. Values in brackets correspond to 
photometric errors in units of valid digits. Upper limits are given at 3cr confidence. 



redwards of the 4000 A break, which allows a reasonable esti- 
mate of M* of a galaxy (e.g., Glazebrook et al.|2004 Ilbert et al. 
[20091 SGL09). 

Systematic uncertainties of up to an average of 0.2 - 0.3 dex 
on Mt are present due to the specific details of the stellar popula- 
tion models and the assumed attenuation/extinction law (|Pozzetti 



et al.'2007 ; Kupcu Yoldag et al. 2007 Kajisawa et al. 2009| Ilbert 



et al. 2010 ). Despite the small sample size of GRB hosts, there 



is evidence for an offset of around 0.2 dex between the presented 
method and the one of SGL09 as shown in Fig|4] In the follow- 
ing the recalculated stellar masses of the long GRB hosts with 
the photometric data compiled in SGL09 are used as a compari- 
son sample. 

The SPS fit returns not only the luminosity and mass of the 
galaxy, but also other physical properties of the host, such as the 
age of the dominant stellar population r, its color-excess E(b-v) 
and the star-formation rate (SFR) derived from the rest-frame 
UV flux. Reported physical host properties are the median of 
the probability distribution of the total grid over all galaxy tem- 
plates at the fixed spectroscopic redshift. Errorbars represent the 
maximum and minimum value of the respective parameter in the 
global ^^-distribution of the multi-dimensional parameter grid. 
Typically, errors are asymmetric and dominated by the uncer- 
tainty of the color excess in the host galaxy, which results in 
logarithmic errors on all galaxy parameters. Absolute magni- 
tudes and masses are compared against the redshift dependent 
galaxy luminosity functions (Willmer et al. 2006; Marchesini 
et al.||2007|l, and stellar mass function ( |Marchesini et al.||2009j 




lof 



9 

j^,jSGL09 



10 



11 



Fig. 4. Comparison of stellar masses from the host sample of 
SGL09. The plot shows host masses which were derived from 



Ilbert et al.|2(JT0l i 



the photometry compiled in SGL09 and following Sect. 4.4.2 
with models from Bruzual & Chariot' (2003 ) against values di- 
rectly taken from SGL09. Error bars are the maximum and mini- 
mum value of the stellar mass i n the g lobal -distribution of 3 x 
10* galaxy templates (see Sect. 4.4.2 1. The solid lines represents 
equality, and the dashed line the median offset. Increasingly 
grey shaded areas show dispersions of 0.2, 0.5 and 1 .0 dex re- 
spectively. The inset shows the distribution of mass differences, 
which has a median of around 0.2 dex. 



5. Results 

5.1. Afterglow extinction and metals-to-dust ratios 

Out of the total of eight afterglows in the sample, the extinc- 
tion properties for five of them (GRBs 070802, 080605, 080607, 
080805, 090926B) are extensively discussed in previous works 



([Kriihler et al |2008[|Enasd6ttir et al.|20"09l|Perley et al.|20I lb 
[Greiner et al.[[20lTf [Zafar et al.[|2011[ l. In these cases, also the 
published afterglow analysis is comparable to the approach of 
this work, and it is therefore not repeated here, with results from 
the literature being summarized in Table [5] For the remaining 
three events (GRBs 070306, 081109, 100621 A), we present ei- 
ther new data and their modehng (GRBs 081109, 100621A), or 
a new analysis (GRB 070306) in the following. 



et al. 



The afterglow of GRB 070306 was discussed in Jaunsen 
( 2008| l. As their analysis is significantly different from 



our approach, we refit the available afterglow data following 
Sect. 4.4.1 The broad-band SED in Fig. [5] is reasonably well 
fit (x^ = 123/108 d.o.f.) with a single power-law continuum 
with spectral index /? = 1.00 ± 0.07, an A™^ = 5.5+|f) mag and 



Nh,x - 2.5_Q 2 X 10 cm at 90% confidence, which imphes a 
metals-to-dust ratio of NH,x/Ay^^ - 4.4 x 10^'cm^^/mag. Given 
the redshift of z = 1.496, and the sparse wavelength coverage 
in the NIR (probing the rest-frame optical red-ward of 400 nm, 
where there is little distinction between local extinction laws), all 
local dust models provide equally good fits to the data of course, 
and within errors compatible values of /3, A™^, and Nh,x- No 
strong statements can be made either with respect to a possible 
presence of a break between NIR and X-ray data. We adopt the 
model with the least number of free parameters (single power- 
law continuum, which also yields the lowest x^), but note that in 
the case of a break between the two wavelength regimes (as seen 
in most early GRB afterglows, Greiner et al.|201l] l the fit is of 
comparable quality (A^^ = 3), and yields an best-fit Ay™ which 
would be significantly lower (A^"^^ = 4.3;^[ q mag), but within 
errors consistent with the single power-law values. 

The afterglow SED for GRB 081109 at z = 0.979 has been 
constructed from GROND and Swift data and is shown in Fig. [6] 
After subtraction of the late host epoch, no residual flux is de- 
tected in the two bluest GROND filters g' and r'. The after- 
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Table 5. Afterglow properties 



Afterglow 


Redshift<*> r'^*' 




4 GRB 


Nh.X 


A-'/d-o.f. 




[s] 

1 An^ -iirv(a) 




[mag] 

c c + 1.2 


10^'cm-- 

o c+O.J 





GRB 070802*'^' 
GRB 080605* '^' 
GRB 080607<''-f* 
GRB 080805* '^' 

GRB 081109 
GRB 090926B"'' 
GRB 100621 A 



2.452 
1.640 
3.036 
1.505 
0.979 
1.24 
0.542 



16.4('" 

20<'=> 
79(g) 
78(h) 

190*" 
110ij> 
63.6« 



0.60 
0.67 
0.96 
0.47 
1.10 
0.73 
0.79 



1 10+""' 
0.67 + 0.01 

0.97 ± 0.05 

|-| 7Q+8:8? 

^■'-'-0.07 
1 9Q+0.11 
A -^^_n in 



■^-1.0 

0.47 ±0.03/ 1.20!" 



1.23^1^/1.19: 



-1+0.09 



2.0+°-' / < 2.9 



-0^0 

3.3±0.3/2.33!;j™ 
1.01!jJ[^/1.53±0.i3 
3.4+°'' 
1.4+ 
3.8+ 



1 01+0.09 

2, 



+0.09 /o 71 + 
!.7+^-^ / 3.8+!;-2 



08 



- ^1 

s;l!:2 
-0.2 



1 0+""' / 1 2+"' 

1 10+"" 

9 9+0.5 

1 62+°'5 



15/14 
428 / 327 
70/39 
23/18 
48/63 
33/31 
138 / 124 



Redshifts froniiJaunsen et al. 


j2008|, 


Eliasdottir et al. 


2009J, 


Fynbo et al. 


2009a|, 


Fynbo et al. 


|2009b 


1 and 


Milvang- Jensen et al. 


j2010| 



2 J ^Kouveliotou et al.|1993) . All GRBs in this work are hence unambigously long events, 
the fitting, and hence has the same error 



long GRBs have Tgo > '. 
'"''/^opt is tied to fix in the 1 

Notes: Afterglow measurements were taken from the reference denoted with the superscript in the first line. When two references are given, we 
quote both val ues for Ay™ and N hj;, but values for /3 andx^ only from the first one for the sake of clarity. 



References: (a)|Barthelmy et al.||2007ll, (b)|Greiner et al.|([20TT]l, (c)|Zafar et aTlpOTi), (d)|Cummings et al .'( 2007), (e) Cummings et al. (20081, (f) 
Perley et al.|(|2011b| l, (g) |Stamatikos et al.||2008^ , (h) |Palmer et ar||2008p , (i) |Markwardt et al.| ( |2008^ ,"(]) |Baumgartner et al., ( ,2009j , (k) ,Ukwatt"a| 
et al.| ( |2010a| > 
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Fig. 5. NIR to X-ray SED of the afterglow of GRB 070306 ob- 
tained approximately 120 ks after the trigger in the observer's 
frame. The dashed line shows the unabsorbed synchrotron con- 
tinuum emission while the best-fit model (including dust and 
metal absorption) is shown by solid lines. Upper limits are 
shown by downward triangles. 



glow is detected in all five redder bands, which implies an ex- 
tremely red color of (R - K)ab S 6 mag and /3ox < 0.44. 
The combined GROND and XRT data are well fit with a sin- 
gle power-law continuum, indicating that both the optical/NIR 
and the X-ray regimes probe the same part of the synchrotron 
spectrum. The obvious curvature in the GROND data is accu- 
ratly described with either of the local dust models, with best- 
fit parameters of A^'^'' = 3.4+[;-4mag, /3 = 1.12 ± 0.02, and 
Nh.x = (1.1 + 0.1) X lO^^cm-- at 90% confidence and a of 
48.1 for 63 d.o.f. in a MW-like extinction law. SMC and LMC 
models yield within errors comparable parameters and provide 
equally good fits to the data (x'& of 48.2 and 48.8, respectively). 

The SED of the af terglow of GRB 100621A at z = 0.542 
( Milvang- Jensen et al.|2010 ) is shown in Fig. [7] Similarly to the 
SED of GRB 081 109 there is strong curvature and obvious red- 
dening in the optical/NIR part of the SED. The inferred ultra-red 
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Fig. 6. S ame as Fig . [sjfor the afterglow of GRB 08 1 1 09 obtained 
approximately 60 ks after the trigger. 



color of (R - K)ab ~ 5.8 mag, and the /3ox value of 0.39 provide 
evidence for strong dust extinction. The best fit is obtained with a 
broken power-law with spectral indicesySi =/?2-0.5 = 0.79+y ],j, 
as well as Ay^^ = 3.8 + 0.2 mag for an LMC-like extinction 



,-2 



law, and Nh.x - (1.62 
(X^ of 138 for 124 d.o.f). Given the rest-frame coverage of 
~ 300 - 1500 nm all local dust models return comparable values 
with A'^™ values of 3 .8 + 0.2 mag for an SMC- and 4.0 ± 0.2 mag 
for a MW-type extinction law. All data bluewards and including 
the r' filter are consistent with this extinction laws, while the 
g'-band photometry is somewhat (2 - 3cr) brighter than the best 
fit predicts. This could indicate a discrepancy between the de- 
tails of the specific dust extinction law and local models similar 
as observed e.g., in GRB s 070802 or 080607 (Eliasdottir et al. 
|2009l|Perley et al.|2011b| l. 

The visual extinctions for GRBs 070306, 081109 and 
100621 A are among the largest ever measured directly along 
GRB sight-lines, and further imply metals-to-dust ratio of 
NH,x/Ay^^ ~ (3 - 5) X 10^'cm-^/mag. Compared to previous 
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Fig. 7. Same as Fig. [s] for the afterglow of GRB 100621 A ob- 
tained approximately 7.6 ks after the trigger 



afterglow measurements, these NH,x/Ay™ values are relatively 
low, and within a factor of 2-3 similar to Nh/Av as observed in 
the SMC, LMC and MW (see also SecL^^. The results of af- 
terglow SED fitting, as well as values taken from the literature, 
are summarized in Table |5] 

5.2. Host properties 

5.2.1 . The host of GRB 070306 

The galaxy hosting the strongly extinguished GRB 070306 at 

( |2008l l. 



Jaunsen et al. 



z - 1.496 was previously discussed in 

In addition to the public VLT imaging data (FORS R, ISAAC 
JsHK,), the host of GRB 070306 was observed with GROND 
(g'r'i'z'JHKs simultaneously), and its SED (Fi g. [8]) is further 
complemented by published u and /-band data ( Jaunsen et al 
2008| . The host is bright (/ = 23.1 mag), mildly rec|]with (R - 



K)ab ~ 1 .5 mag and shows evidence of a 4000 A break. The data 
are well fit with a non-extinguished (Ay < 0.6rlhost template, 
and yield an absolute AB magnitude of Mb - -22.4 ±0.1 mag, 
which, at z ~ 1.5 corresponds to ~ 1.7 L*. The stellar mass of 
log(M»[MQ]) - 10.4 + 0.2 puts the galaxy among the most mas- 
sive hosts compared to the sample of SGL09. The star-formation 
rate estimate from the rest-frame UV flux is 13^1^ A^o/yr, which 
gives a specific star formation rate (SSFR = SFR/M,) of ~ 0.5 
Gyr-', or growth timescale (i.e., 1/SSFR) of 2 Gyr. The SFR 
is in reasonable agreement with the one derived from the [Oil] 
emission line ( [Jaunsen et al. 2008). We note, that the ISAAC 
//-band host image was obtained only 2.5 days after the GRB, 
and hence very likely contains a significant fraction of afterglow 
light, explaining the blue H - K color. The physical properties 
of the host, however, are comparable if fit with or without the 
ISAAC H filter. 



^ We compare GRB host colors against the median (R 
0.8 mag color from the SGL09 sample 

^ For galaxies, we use Ay, or Eg y to indicate measured attenuation 
and effective reddening, since these quantities depend for example on 
the topology of the ISM and galaxy geometry i jPierini et al.|2004) , and 
on galaxy scales are different from the corresponding values of a given 
extinction law. 



5.2.2. The host of GRB 070802 

The host of GRB 070802 at z = 2.452 was discovered in deep 
FORS R and ISAAC K band imaging ( [Ehasdottir eTaL][2009) . 
The afterglow SED is characterized by a prominent 2175 A dust 
feature, and significant dust in the range of Ay ~ 1 mag ( Kriihler 
et al.|2008|[Elfasd6ttir et al.|2009| ). To construct the optical/NIR 
host SED, the R ~ 25.2 mag host is further observed with 
EFOSC/NTT in / and HAWK-IA'LT in the 7-band. The galaxy is 
moderately red ((R-K)ab ~ 18 mag) and its SED (Fig.[8]l shows 
a 4000 A break, but the age of the dominant stellar population is 
not well constrained by the available data (< 1 Gyr). There is no 
strong evidence for internal reddening, and the best-fit absolute 
magnitude is Mg - -21. 4 ±0.2 mag, which is ~ 0.6 L* at z ~ 2.5 
with log(M,[M0]) = 9.7+°-2. Using the rest-fr ame UV flux de- 
rived from the galaxy model fitting, an estimate for the SFR of 
10!^" Mo/yr, and the SSFR of ~ 2 Gyr-' can be derived. 



5.2.3. The host of GRB 080605 

The host of GRB 080605 at z = 1.640 was discovered in late 
GROND follow-up observations of the burst field 22 days af- 
ter the GRB trigger. The afterglow SED shows evidence for 
a 2175 A dust feature, and significant Ay in the range of ~ 
0.5 - 1.3 mag ( [Greiner et al.|2011[|Zafar et al.|[20n] i, and will 
be further discussed in Nicuesa et al. (in preparation). The host 
is bright (r' ~ 22.8 mag), and blue with a flat g' - z' color, and 
(R - K)ab ~ 0.5 mag as estimated from the best fit galaxy model 
(Fig.[8}. The SED fit further yields Mb = -22.6±0.2 mag, which 
is ~ 2.1L* at z ~ 1.5, and log(M,[MQ]) = 9.6+°-^. The domi- 
nant stellar population of the host is young (r = 0.06;'^y Q3 Gyr), 
and there is no evidence for reddening at the 2cr level (E'^_y ^ 
0.22 mag). The host galaxy is vigorously star-forming with a 
SFR of 40+^° Mo/yr and a SSFR of ~ 10 Gyr-' . 



5.2.4. The host of GRB 080607 

The afterglow of GRB 080607 is heavily reddened (Ay ~ 
3 mag), has a modest 2175 A dust feature and is characterized by 
a strong neutral hydrogen absorber with roughly solar metallicity 
and molecular gas (Prochaska et al.|2009 ; Perley et al. 2011b} . 



Data from Chen et al. ( 201 1 1 show aR ~ 27 mag, very red host 
with a synthetic color of (R - K)ab ~ 3 mag (Fig. [S]). The host 
is well described with an extinguished galaxy template (Ay ~ 
1 - 2 mag), and physical parameters of Mb - -21.1 ± 0.3 mag. 



which is ~ 0.3 L* at z ~ 3, log(M,[M0]) = 9.9!°-^, an extinc- 
tion corrected SFR of 40+^^ MJyi, and SSFR of ~ 8 Gyr-'. 
These values are in good agreement with previously published 
ones ( |Chen et al.|2011| l. 

5.2.5. The host of GRB 080805 

GRB 080805 had a very red afterglow, where both an SED 
and spectral analysis showed a large dust column (Ay ~ 1 .0 - 



1.5 mag), and evidence for a 2175 A dust feature (Fynbo et al. 
2009a, Greiner et al. 2011; Zaf ar et aLllSoTT] !. Its host is dis- 
covered in late EFOSC/HAWK-I imaging in five filters (VRiJK, 
see Fig. [8]l, is relatively bright (R ~ 25.5 mag) and red ((R - 
K)ab - 2.5 mag), with best-fit physical parameters of Mb = 
-20.4±0.2mag(~ 0.3 L* atz ~ 1.5) and log(MJM0]) = 9.7+j^j. 
The remaining galaxy properties are not well constrained by 
the available data, yielding a limit on the galaxy reddening 
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Fig. 8. SEDs of the host of in this sample and the best-fit galaxy model (soUd Une) in the observer's frame. Filled black circles 
represent photometric measurements, while downward triangles denote 3cr upper limits. 
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Table 6. Host galaxy properties 



Host 


^spec 


■Sphot 


R.A.*''' 


Decl.<°> 


Offset**' 


Probability*''' 


Instruments*''' 


Filter 








J2000 


J2000 


arc sec 








GRB 070306 


1.496 


1 64+"^^ 


09:52:23.31 


+ 10:28:55.4 


< 0'.'2 


0.002 


FGl 


g'r'Ri'z'JaJHK^" 


GRB 070802 


2.452 


2 3+'-' 


02:27:35.72 


-55:31:39.0 


^.'15 ±a'08 


0.011 


EFHl 


RIJK 


GRB 080605 


1.640 


irS:^ 

I. / „ 2 


17:28:30.05 


+04:00:56.0 


a'23 ±a'ii 


0.002 


G 


g'r'i'z'JH 


GRB 080805 


1.505 


1 7+0.6 
^•'-0.7 


20:56:53.43 


-62:26:39.3 


<0'.'3 


0.030 


EH 


VRIJK 


GRB 081109 


0.979 


0.84!«;i« 


22:03:09.63 


-54:42:39.9 


< a'2 


0.002 


EFGHU 


Ug'Vr'i'Iz'YJaJHK 


GRB 090926B 


1.24 


J 43+0.62 


03:05:13.91 


-39:00:22.6 


<a'6 


0.018 


EFGS 


Ug'r'i'z'JHK 


GRB 100621 A 


0.542 


50+^-^^ 


21:01:13.08 


-51:06:22.2 


a'i2±a'08 


0.0006 


EGHU 


wlmlw luUgY i'z' YJHK 



Notes: Redshifts from |Jaunsen et al.H2008) , |Elfasd6ttir et al.H2009^ , |Fynbo et al.| (|2009ai, Fynbo et al. (2009b) and ' Milvang- Jensen et al.| ( |2010| l. 
The host of GRB 080607 is not shown in this table. All measurements were directly taken from Chen et al. (2010 20Tr|i 

Host position derived after tying the astrometric solution to the USNO-Bl catalog (Monet et al. 2003 1. Typical absolute uncertainties are a; 0'.'3. 

Relative offset calculated by registering the host images against astrometric templates derived from afterglow images with a typical precision of 
40 mas rms 

Estimated chance coincidence probability following Bloom et al. j2002| and Perley et al. ( 2009 1 

G is GROND at the 2.2 m MPG/ESO telescope, E/S are EFOSC/SOFl at the NTT, and F/HTare HAWK-1, FORSl/2 and ISAAC at the VLT, 
and U UVOT onboard Swift respectively. 



Further complemented by the ii and / band magnitudes in Cool et al. 1 2007 1 and Jaunsen et al. 1 2008 1 



Table 7. SPS Host galaxy fitting results 



Host 


Redshift 


Mb 


Age 


F' 


log(M.) 


log SFR 


log SSFR 


L 








magAB 


Gyr 


mag 


Mo 


Mo/yr 


yr"i 







GRB 070306 1.496 -22.4 ±0.1 1.6+1;; < 0.16 (2a-) 10.39!^; |^ 1.1+""^ -'^■'i a\ 1-7 + 0.2 11.8/12 
GRB 070802 2.452 -21.4 + 0.2 0.38™ < 0.42 (2cr) 9.7+^-^ 1.0+^"° -8.7+™ 0.6 + 0.2 2.3/4 



GRB 080605 1.640 -22.6 + 0.2 0-0(>_om < 0.22 (2o-) 9.6!^"^ l.e^^^^ -^■'^ ol 2.1 + 0.4 3.3/6 

GRB 080607 3.036 -21.1 + 0.1 0.16!!!" 0.35^025 9 9+0:4 j g+M -8.2+i;:^ 0.3 + 0.1 3.9/5 

GRB 080805 1.505 -20.4 + 0.2 0.51!||^ < 0.65 (2o-) 9.7+°| 0.8™"^ "^-^-os 0.3 + 0.1 0.3/5 

GRB 081109 0.979 -21.27 + 0.09 0.24!|jJ 0.24!"!!« 9.82+*™ 1.5+^"^ -8.4!|^ 0.9 + 0.1 11.2/14 

GRB090926B 1.24 -21.5 + 0.1 0.14!^°* 0.35!™ 10-l!of l-9!„, -8.1™ 0.9 + 0.1 7.3/9 



GRB 100621 A 0.542 -20.68 + 0.08 0.05!""' 0.14!"";^ 8.98!"+ l-13!n'n -7.9!"-^ 0.6 + 0.1 16.5/14 



-0.5 -Q.5 
+0.14 1 1 T+0.15 



E'g y < 0.7 mag, an age of the stellar population of 0.53!° 42 Gyr^ 
a SFR of 6!f Mo/yr, and SSFR of ~ 1 Gyr ' . 

At a distance of 2.5" north-east of the afterglow/host posi- 
tion, there is another R ~ 25 mag, and even redder ((R - K)ab ~ 
4 mag) galaxy, a plausible candidate for the strong Mg II absorb- 
ing system at z = 1 .20 reported in Fynbo et al. ( 2009a 1. 



5.2.6. The host of GRB 081 1 09 

GRB 081109 is the only burst in the sample where no spectro- 
scopic redshift was available in the literature. However, the host 
is bright (r' ~ 22.7 mag), moderately red {{R - K)ab -1-6 mag) 
and our host spectrum (see Fig.|9|l reveals a single emission line 
above a well-detected continuum. Within the wavelength cov- 
erage of the spectrum (~ 370 - 950 nm), this emission line is 
interpreted as [Oil] [/13727] at z = 0.9787 ± 0.0005. If it were 
any of the other prominent nebular lines (H/3, [OIII], Ha), we 
would have expected to detect [Oil] [A3127] in the spectrum as 
well. At this redshift, there is further spectroscopic evidence for 
a Balmer break, and a tentative absorption of the Call HK dou- 
blet (see Fig. |9]). 

The host SED is shown in Fig. [8] and well fit with a young 
(t = 0.1 8!{},',^ Gyr) and reddened (A'y = 1.0 ± 0.2 mag) stellar 
population. The absolute magnitude of Mb = -21.27 ±0.09 mag 
corresponds to ~ 0.9 L* at z ~ 1 . The stellar mass obtained from 
the SPS fit is log(M.[M0]) = 9.80!;}j;^, with a SFR of 33![^ 
Mg/yr, which together yields a SSFR of ~ 5 Gyr The emis- 
sion line flux of [OII][/13727] is (1.8 + 0.2) x lO-'^'erg/cmVs, 
which includes a systematic error contribution from the abso- 




10 



600 700 800 

Wavelength [nm] 



Fig. 9. Wavelength and flux-calibrated FORS2 300V spectrum 
of the host of GRB 08 1 109 in black. The thin grey line shows the 
sky spectrum. The red line is the best-fit galaxy model obtained 
from the available photometry and red points are the photometric 
measurements. 



lute flux normalization. This correspo nds to a d ust un-corrected 
star-formation rate of 13 ± 4 Mg/yr ( Kennicutt^l998 ), or 48!}^ 
Mo/yr when using the upper A'y measurement, which is in good 
agreement with the SFR from the UV continuum. 
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5.2.7. The host of GRB 090926B 



Promptly after the trigger, Fynbo et al. ( 2009b| l reported the de- 
tection of the host galaxy of GRB 090926B based on a promi- 
nent [Oil] emission line in a VLT/FORS spectrum. The host 
is also imaged in later GROND and NTT observations with a 
brightness of r' ~ 23 .0 mag, and a mildly red color ((R-K)ab ~ 
1 .5 mag) as shown in Fig. |8] The host SED is well fit with an ex- 
tinguished host model (A'^ = 1.4+g2mag) withMfi = -21.5+0.2 
(~ 0.9 L* at z ~ 1.3) and log(MJMa]) = 10.1!°'^. The SFR 
based on the extinction corrected UV flux is SO;^^^'' Mg/yr which 
is among the highest ever measured from optical data, but not 
well constrained given the uncertainty on the dust extinction 
properties. The constraints on the age of the stellar population 
and SSFR are weak, with values of t = 0.14+°gg Gyr and the 
SSFR of ~ 7 Gyr ', respectively. 



5.2.8. The host of GRB 100621 A 

The host of GRB 100621A was reported very early after the 
trigger as a DSS2 source, providing a constant contribution to 



the afterglow in the UV/blue light curve (Updike et al. 2010 
lOates & Ukwattal[20T0| ). In fact, the redshift of z = 0.542 of 
GRB 100621 A is based on emission lines from a bright host 
( [Mil vang- Jensen et al.|[20T0| . The SED of the r' ~ 21.5 mag 
host is well-sampled from the UV to the NIR and shown in 
Fig. |8] In strong contrast to the extremely red afterglow {R - 
K)ab ~ 5.8 mag, the host is blue with an inferred color of 
(R - K)ab ~ 0.3 mag and {uvw2 - K)ab ~ 0.9 mag. The 
SPS host fit returns an intrinsic extinction of Ay = 0.6^q 2 mag 
for a very young stellar population of age t - 0.05 ^q J^ Gyr 
This host has the lowest stellar mass in the presented sample 
with log(M.[MQ]) = 8.98+;j[4, and an absolute magnitude of 
Mb = -20.68 + 0.08 mag, which is ~ 0.6 L* at z ~ 0.5. The SFR 
and SSFR are 13;^^ MJyr and ~ 14 Gyr"', respectively. 

6. Discussion 

6.1. Dust reddening in GRB afterglow SEDs 

The visual extinction measured from X-ray-to-NIR SED-fitting 
towards GRBs 070306 {A'^^^ = 5.5![;2mag), 081109 (A^^^ = 



3.4+;;-;^mag) and 100621 A (A^'^'^ = 3.8 ± 0.2 mag) are among 
the largest ever derived from optical/NIR data for GRB after- 
glows (e.g., Savaglio et al. 2003 [ Kann et al. 2006[ Greiner 
|et al.|[20TT] l. They clearly show that a large dust column can 
be the dominating feature in a GRB afterglow SED. The dust 
properties inferred from afterglow measurements are well rep- 
resented with local models in the rest-frame 300-1100 nm, and 



at the resolution of broad-band imaging (see also e.g., Fynbo 



et al. 


|2001 Watson etal. 


2006 Kann et al.|2006 Schady et al. 


2007 


IStarhng et al. 


2007 


, while noteworthy exceptions exist in 


the literature (e.g.. 


SavagUo & Fall 2004 Perley et al. |2010b 



Clemens et al. 201 l| l. The good fit provided by local dust ex- 



tinction laws further suggests an abundance of small dust grains 
comparable to the MW/LMC or SMC. There is hence no di- 
rect evidence that the dust towards these GRBs through their 
hosts is different than observed along local sight-lines. A differ- 
ent dust grain size distribution would have been expected if the 
dust were located in the immediate vicinity {R < 10'^ cm) of 
the GRB and shaped through its intense radiation, i.e., through 
dust destruction (Waxman & Draine 2000 ; Fruchter et al.|20()T| 



Nh/Av ~ 2 X lO^'cm 2/mag ( |Predehl & Schmitt|[T995| l. For 
un-extinguished GRB sight-lines this ratio is generally found to 
be a factor of 10-100 times larger than those of the Magellanic 
Clouds or Milky Way (e.g., |Galama & Wije rs 2001 ; Strat ta et al.| 
2004) . This is suggestive of a dependence of the metals-to-dust 
ratio on the amount of visual extinction. We will return to this 
issue in Sect. 16.41 



6.2. The hosts of the dustiest afterglows 

The general properties of the selected GRB host galaxies are di- 
verse. They have (R - K)ab colors ranging from flat and blue 
(R - K)ab ~ mag to extremely red {R - K)ab ~ 3 mag with 
an average color of {(R - K)^^^} - 1.6 mag, and host extinction 
values between Ay ~ Omag and Ay ~ 2 mag. Also, their stel- 
lar mass and absolute magnitude distributions are broad, with 
values between log(MJMQ]) = 9.0 to log(MJMo]) = 10.4 
«log(M.[M0])) = 9.8 + 0.4) and Mb between -20.3 mag and 
-22.6 mag {{Mb) - -21.3 ± 0.6 mag). These absolute bright- 
nesses are in a range between several tenths to few L* ((L) = 
0.9L*) as compared to the general field galaxy population at the 
same redshift. The average SFR and SSFR are about 30 M^/yr 
and (log SSFR [yr]) ~ -8.3, respectively. The average growth 
time is ~ 0.2 Gyr, which illustrates that not only optically- 
selected hosts, but also those of highly-reddened afterglows are 
very efficient in producing stars. A rough estimate on the metal- 
licity of the hosts can be obtained if these GRB hosts follow 
the fundamental plane as defined from nearby SDSS galaxies 
(Mannucci et al. 2010a). With a given stellar mass and SFR, the 



host galaxies in this sample are expected to have metallicities in 
a range between 12 + log(0/H) ~ 8.2 and 12 + log(0/H) ~ 8.9, 
with an average of 12 + log(0/H) ~ 8.6. We caution that the 
SFRs were derived using the rest-frame UV flux, which is quite 
sensitive to the dust extinction properties. 

Although not well-constrained in all cases, the average 
luminosity- weighted effective reddening inferred from host pho- 
tometry is typically smaller or equal to that measured from af- 
terglow observations. This is not a particularly surprising result, 
given that the sample selection was based on high visual extinc- 
tions of the afterglow SEDs in the first place. It directly indicates 
some variation in the dust distribution of the hosts, which again 
is not a surprising result, given the differences in extinction prop- 
erties along different sight-lines through the difiiise ISM to Giant 
Molecular Clouds in the Local Group (e.g., Gordon et al.|200J 



[Draine & Hao||2002 i. In addition, the metals-to-dust ratios for 
these afterglows are only a few times the Galactic value of 



Fitzpat rick|2004| l, and the geometrical differences between a sin 
gle sight-line and an extended distribution of star-light and dust 
(e.g., |Gordon et al.|1997[[Silva et al.|1998) . 

One intriguing case is the host of GRB 100621A. Although 
having one of the most extinguished afterglows ever detected 
(even in the presented sample), its host shows very blue colors, 
and is one of the youngest and the least massive galaxy in this 
work. This particular example provides evidence for a patchy 
dust component where the geometry of the dust distribution and 
not the properties of the host galaxy makes the single GRB sight- 
line dust-enriched. 



6.3. Comparison to previous GRB host samples 

One key result of this study is the success rate of the discovery 
of the selected hosts. Out of eight hosts, which were selected 
given their afterglow properties (hence a selection independent 
on host properties, in particular galaxy brightness), all are lu- 
minous enough to detect in optical ground-based imaging. This 
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fraction is significantly larger than expected from a host sam- 
ple based on XRT detections ( Fynbo et al.|2008 Malesani et al. 
|2009 ). The effect of an increased detection rate is even stronger 
in the NIR: Seven out of eight are detectable in the /T-band, while 
this fraction is only ~35% for the general host population of 
Swift/GRBs ( [Malesani et al.|2009| l. Partially, this is the result of 
the lower average redshift of the selected hosts ({zav) - 1-5) as 
compared to all Swift GRBs with (zswifi) ~ 1-9. 

The lower redshift is however not the only reason for the high 
detection rate. The selected hosts are on average redder and, as 
shown in Fig. [TO] have typically higher luminosities and stel- 
lar masses than the (sensitivity-limited) SGL09 sample which 
has {{R - K)ab} = 0.8 mag, {Mb) = -19.6 ± 1.5 mag and 
<logM,[Mo]) = 9.1 + 0.6. A two-sample K.-S. test returns p- 
values of 0.002 for the stellar mass, and 0.006 for the absolute 
magnitude distributions respectively, which is tentative evidence 
that both distributions are not drawn from the same parent sam- 
ple. However, given the small sample size of only eight high-Av 
events, larger samples are required to statistically establish the 
existence of a difference at higher significance. Of course, both 
distributions are drawn from the same physical parent sample 
(GRB hosts), indicating that the different selection criteria probe 
different host properties. 

A possible explanation of the different host properties would 
be the now on-average higher redshift relative to the SGL09 sam- 
ple, where star-formation was driven by more massive galax- 
ies as compared to the more nearby Universe (e.g., |Cowie et al. 
1996 Hopkins 2004 1. To test this hypothesis, we selected a sub- 
sample from SGL09 with a median redshift comparable to the 
hosts in this work. This essentially removes all z < 1 SGL09 
hosts and leaves only 13 events for comparison (see histograms 
in Fig 
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Despite the small number statistics, the and Mg 
values are again placed at the high-mass and high-luminosity end 
of their respective distribution, and a K.-S. test is also marginally 
suggestive of a difference (/^-values of 0.001 and 0.034 for the 
masses and absolute magnitudes). 

We conclude that by selecting extinguished afterglows 
we are very likely probing a more luminous, massive and 
chemically-evolved population of GRB hosts. 

As it is clear that these were largely missing from previous 
samples due to their poor localizations, there is a selection bias 
and the host population is missing most of its massive, evolved 
and metal-rich members. As a direct consequence, GRB hosts 
trace the global SFR closer than indicated in studies which are 
based on host samples of optically selected GRB afterglows, and 
the apparant deficiency of high-mass host galaxies is at least par- 
tially a selection effect. 



Similar conclusions apply for all galaxies hosting afterglows 
that show a significant 2175 A dust feature in their SED. Four 
out of five cuiTently known afterglows are within the presented 
sample, which argues for a direct connection between large dust 
columns and the presence of the UV bump (see e.g., | Greiner| 
et al. 201 1[ Zafar et al. 201 1\ . Their, on average, more mas- 
sive and luminous hosts suggest a qualitative relation between 
the stellar mass of a galaxy and the presence of a 2175 A fea- 
ture, where the latter is only present in fairly massive and 
metal-enriched galaxies (see also e.g., Noll & Pierini| 2005 1. 
Conversely, a strong 2175 A feature in an afterglow SED is also 
very likely a good proxy for the stellar mass and luminosity of 
the GRB host. 
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Fig. 10. Distribution of stellar masses and luminosities of the 
hosts of highly extinguished afterglows (blue) and the host sam- 
ple from SGL09 (grey). 
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Fig. 11. Distribution of stellar masses and luminosities of the 
hosts of highly extinguished afterglows (blue) and a subsample 
of SGL09 (gi-ey) with <z) -1.5. 



6.4. Metals-to-dust ratios in context 

The ratio between the line-of-sight extinction and the total metal 
column for GRB afterglows has been investigated in a num- 
ber of papers (e.g., |Galama & Wijers[ 2001 ; S tratta et al 



ISavaglio & Fall||2d04'; 'Kann et al.||2b06( |Schady et al. 



2004 



2007 



2010; Greiner et al. 2011; Zafar et al.||2011 l where ratios typ 



ically much higher than the ones observed in the Local Group 
were derived. Measurements fbr different Galactic s ight-lines 
(e.g.. 



Predehl & Schmitt 



almost universal value of 



1995 



Giiver & Ozel 



|2009| l show an 
■ around Nh/Av * 2 x lO^'cm^^mag, 
while the matter probed by afterglows can yield metals-to-dust 
ratios up to and sometimes even above 100-times higher (e.g., 
[Watson et al.|2007|[Rairet al.|2010| ). 
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Fig. 12. Metals-to-dust ratio versus sight-line extinction for GRB 
afterglows. The horizontal dashed lines marks the selection cri- 
terion for GRBs to enter this sample. Vertical dashed lines il- 
lustrate metals-to-dust ratios of 1, 10 or 100 times the Galactic 
value. Solid lines show the toy model of two physically inde- 
pendent absorbers, where one is fully devoid of dust with a H- 
equivalent metal column of Nh.x - 10^'^cm"^ and represented 
by the dotted line, while the other is neutral and has metals-to- 
dust ratios of 1, 2 or 3 times the Local Group value. Two in- 
dividual cases (GRBs 080605 and 080607) illustrate the scatter 
between the analysis of different data sets. 



6.4.1. An anti-correlation between metals-to-dust ratio and 
sight-line extinction 



Figure 12 shows the Nh,xI^'v^^ ratio for a large number of af- 
terglows and illustrates its dependence on the sight-line dust ex- 
tinction. With the afterglows in this work, there is now for the 
first time reasonable coverage in the A^*^^ ~ 1 - 5 mag range. 
Intriguingly, the metals-to-dust ratio is strongly anti-cor related 
with A^*^^, confirming the tentative trend reported by Perley 
et al. (2009 1. A Spearman rank-order correlation analysis tor tlie 
combined sample in Figure 12 returns a correlation coefficient 
p — -0.63, with a two-tailed ;?-value of 3 x 10"^, in strong con- 
trast to a constant, universal NH,x/Ay™ ratio. 

There are two straight forward ways to reconcile this result: 
One is a dependence of the metals-to-dust ratio on the specific 
environment such that evolved and dust-enriched hosts are more 
efficient in forming dust out of their metals (and in fact we do 
observe on average larger stellar masses for the hosts of high 
Ay^^ afterglows). The other is the presence of two physically 
independent absorbers, where the first dominates the total metal 
column, the other the visual extinction measurements. This triv- 
ially produces a non-correlation between Nh.x and Ay^^, and 
consequently an Nn.x/Ay^^ to A™^ anti-correlation. 

Some outliers of the metals-to-dust anti-correlation (Fig. [T2] l 
might be explained with difficulties of measuring the respective 
physical parameters. This is also illustrated by the example of 
two individual events (GRBs 080605 and 080607) where differ- 
ent values have been published in the literature. Assumptions on 
the continuum emission and the extinction law, notably the total- 
to-selective reddening Ry, but also its parametrization can affect 
the Ay^^ measurement. In addition, there is the possibility that 
the cooling break is located close to or within the range of the 
UV/optical/NlR measurements. In a standard analysis, the intro- 



Fig. 13. Total metals-to-dust ratios for GRB afterglows versus 
stellar mass of their host galaxies. Black data are galaxies host- 
ing a highly-extinguished afterglow using the first values of 



Table 5 where the three hosts with the highest Ay are la- 
beled, i he shaded area indicates the probability distribution for 
optically-selected GRBs, re presented by a log-n ormal distribu- 



tion of NH,x/Ay^^ based on Schady et al. 



(2010 1 and a Gaussian 



distribution in log(MJM0]) based on SULuy. 



duced curvature caused by the spectral break is then interpreted 
as an increased dust column ( Kriihler et al. 201 Ij ). The Nh.x mea- 
surements are prone to errors as well: spectral variation intrinsic 
to the afterglow can lead to wrong estimates on the soft X-ray 
absorption (e.g., Butler & Kocevski|2007[ l. 



6.4.2. Metals-to-dust ratio compared to host mass 



As shown in Sect. 6.2 the hosts of dusty afterglows are on av- 



erage more massive and luminous than their non-extinguished 
counterparts, but there is a broad range of galaxy properties and 
the only common feature between all afterglows/hosts in this 
work is hence the dusty line of sight. In particular, if the envi- 
ronment were responsible for the observed Nh.xIA'^^ to visual 
extinction anti-correlation, we would expect the metals-to-dust 
ratio for GRB 100621 A to be comparable to the bulk of optically 
bright afterglows. It is, however, one of the lowest ever observed 
for GRB afterglows and a factor of 5 lower than the median for 
afterglows with hosts of similar mass (see Fig.[T3]l. Although we 
note that number counts are still too low to derive strong con- 
straints with high statistical significance, this suggests that the 
specific host environment is not responsible for the observed de- 
pendence on the metals-to-dust ratio to the visual extinction. 

6.4.3. A system of two absorbers 

In the second scenario, two to first order physically independent 
(one neutral, host-galaxy related, one ionized, circumburst spe- 
cific) columns of material contribute to the observed absorption, 
denoted as A^//,neutrab A^ff.ion in the following. Nh,x measures the 
sum of both, whereas the A^^^ column would only be associated 
with the A^/f,neutrai absorbcr. Such a circumburst environment is 
not unexpected: The intense afterglow radiation should not only 
photo-ionize the vicinity of the burst, but also destroy the asso- 
ciated dust to large amounts ( Waxman & Draine 2000] |Draine & 
|Hao|2002|[Perna et al.|2003] l, albeit with different effective radii. 
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For A^"^^ ~ 0.1 mag sight-lines, A^n.neutrai A^ff.ion directly 
results in a large metals-to-dust ratio. With an increasing Ay™ 
column of around 1 mag, the A^//,neutrai absorber contributes sig- 
nificantly to the total metal column (A^H.neutrai ~ A^H.ion): this is 



illustrated in Fig. 12 by the solid lines, where the NH,x/Ay ra 



tio asymptotically reaches 1, 2 or 3 times the Local Group value. 
Even larger Ay™ columns than present in this work can test 
this hypothesis. For an Ay™ > 10 mag sight-line, for example, 
A^H.neutrai IS much larger than Nn.ion and the expected metals-to- 
dust ratio would be comparable to its intrinsic value, and to the 
Local Group value (if universal). 

For the bulk of standard, un- or mildly extinguished after- 
glows, the large column of ionized metals with an equivalent 
NH,ion ~ 10^'"^^ cm"^ in the circumburst material dominates the 
total Nh,x measurement, whereas the visual extinction is very 
likely caused by dust further out, either in the diffuse ISM, or 
localized in interstellar clouds (see also Sect. |6.5| l. Hence, for 
optically bright GRBs the total metal absorption as probed by 
the soft X-ray absorption is neither a good measure nor a di- 
rect tracer of the dust extinction along the line of sight. Or, in 
other words, an Nh,x column even as large as 10^^ cm"^ does 
not necessarily imply a significant visual extinction (see also, 
e.g., [Galama & Wijers|[2d&T] |Kann et al.|[2006l [Schady et al. 
|2007| i; 

Also, the shape of the extinction law, which is typically con- 
sistent with sight-lines through galaxies of the Local Group ar- 
gues against an association of the absorbing dust with the im- 
mediate vicinity of the burst. The steepness of the UV rise of 
the extinction law as probed by GRB afterglows is comparable 
to the one of the LMC ( Schady et al.|20irb| , generally sugges- 
tive of an abundance of small dust grains. These grains, however, 
are expected to be destroyed first and the grain-size distribution 
function would be skewed to larger grains. This would result in 
flat, or grey extinction properties, in contrast to the ultra-red af- 
terglows of GRB 081109, or GRB 100621 A, for example. 

We thus conclude that the anti-correlation between metals- 
to-dust ratio and sight-line extinction indicates the presence of 
two absorbing systems, which are to first order physically inde- 
pendent. One of them is dusty, the other ionized and dust-free, 
where the former is probed by the optical/NlR data and the latter 
typically dominates the Nh.x measurement. 



6.5. Location and geometry of the absorbing dust column 

A natural question about the nature of the absorbing dust, gas 
and metal columns detected in the afterglow SEDs and spectra 
is their locations, and if they are directly related to the burst en- 
vironment. A number of previous studies have already revealed 
some clues about the geometry of the absorbing matter, which 
have been derived quite exclusively from sight-lines with low to- 
tal dust content: The distance of the cold-neutral material, linked 
to the DLA and the low-ionization metal absorption lines has 
been constrained to few hundreds of pc to few kpc ( [Prochaska 
|et al.||2006t [Vreeswijk et al.||2007| ). In contrast, the high metal 
column densities as derived from soft X-ray absorption were as- 
sociated with a fully ionized circumburst medium up to few to 
several tens of pc (Watson et al. 2007; Schady et al. 201 la). For 
bursts with largely unextinguished afterglows, most of the met- 
als along the sight-line are typically in a highly ionized state and 
only <10% of the absorbing gas is neutral (Schady e t al.|2011a) l. 
In addition, there is no statistically significant correlation be- 
tween the soft X-ray absorption and the dust column as inferred 
from optical/NIR data nor the metallicity of the neutral material 



(e.g., [Schady et al.||2007| |Zafar et al.|[20n] ), nor the darkness 
of a GRB fCampana et al. 2010[ l. There is however a trend of 
higher visual extinctions with larger neutral metal columns, an 
anti-correlation between gas-to-dust ratio and metallicity (Zafar] 
et al.|20rij ), and dark bursts have stronger neutral metal absorp- 
tion lines in their optical spectrum ( Christensen et al.|201 1[ ). 

The dustiest afterglows and their hosts add two further hints. 
Firstly, there is the previously discussed anti-correlation between 
metals-to-dust ratio and sight-line extinction, albeit with a large 
scatter for individual events. And secondly, their hosts are on av- 
erage redder, more luminous and massive, and supposedly also 
more evolved, dust- and metal-rich than t heir low-A^*^^ co unter- 
parts (e.g., |Kobulnicky & Kewley||2004{ [Savagho et aL][2005 



There is hence a relation between the dust along the sight-Une 
towards the GRB and physical host properties. Such a relation is 
expected if the dust probed by the afterglow is located at large 
enough distances to be fairly representative of the size of the 
host galaxy, its global dust enrichment and chemical state. In 
contrast, dust in the surrounding environment of a GRB, even if 
present and able to survive the intense afterglow or progenitor 
radiation, would rather be related to the very specific details of 
its circumburst environment or the GRB's birth cloud. 

An overlap with optically selected hosts clearly exist: very 
young, blue and low-mass galaxies possess sight-lines through 
dusty regions as demonstrated through the afterglows and hosts 
of, e.g., GRBs 080605 and 100621A. For this kind of events 
the data are clearly not consistent with a uniform dust-shield. In 
a few cases the observations rather indicate a patchy dust dis- 
tribution where the dust is located in clumps of small enough 
covering factor and large enough extinction to be negligible in 
the integrated host light distribution. The variations on the vi- 
sual extinction in some cases are therefore a geometrical effect 
(see also, e.g., |Berger et al.|2003[|Perley et al.|2009l ). 

This raises questions about the validity of star-formation es- 
timates obtained for young and blue galaxies. Afterglow obser- 
vations show that there is extinguished star-formation ongoing 
even in apparently unextinguished galaxies. Here, either the cov- 
ering factor of the dusty clump is just too small to remove signif- 
icantly from the host light, or the clump completely absorbs all 
UV light from the star-forming region, thus having a negligible 
effect on the host galaxy colors. In the latter case, the UV-derived 
SFRs would strictly represent lower limits on the ongoing star- 
formation in these galaxies. Far-infrared and sub-mm observa- 
tions with Herschel and ALMA, respectively, would enable to 
measure the galaxies SFR, dust mass and temperature and would 
help to clarify this issue. 

A coherent picture of highly extinguished afterglows in com- 
bination with their diverse, but on average redder and more mas- 
sive hosts could be obtained within a complex dust distribution 
made out of several constituents related to extinction in the dif- 
fuse ISM, extended interstellar clouds, or localized in fairly com- 
pact and dense regions such as giant molecular clouds. 

The dust extinction is hence very likely not directly associ- 
ated with the GRB environment but plausibly with the neutral 
absorber at distances of a few hundred pc to one kpc. We stress, 
however, that the effective radii of dust destruction and photo- 
ionization will shape the detailed gas-to-dust and metals-to-dust 
ratios, adding further complexity to the absorbing system(s) in 
front of GRB afterglows. Furthermore, the dust distribution in 
high-z galaxies could be even more complex due to dusty galaxy 
outflows, and reflect the absence of a uniform chemical enrich- 
ment on scales up to one kpc (e.g., |Noll et al.|20^9| . 
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7. Conclusions 

The afterglows of GRBs 081109 (A^^^ = 3.4+°-4mag) and 
100621A (A'^^^ = 3.8 ± 0.2 mag) join the growing sample of 
highly extinguished events. Their continuum emission is well- 
constrained by the combination of X-ray and NIR data, and the 
optical observations provide a detailed measurement of the dust 
properties along the sight-lines. While some diversity in their 
extinction properties, particularly dust abundance, clearly ex- 
ists, GRBs 081109 and 100621 A provide compelling evidence 
that a highly obscured afterglow is also a highly reddened one, 
and that extinction laws derived from local sight-lines accurately 
estimate the dust properties towards even highly extinguished 
GRBs. 

Having a large enough sample of coeval afterglows with 
multi-wavelength data would ideally enable to advance from sin- 
gle sight-line, pencil beam investigations to a statistically sym- 
metric geometry where each GRB afterglow represents a differ- 
ent sight-line through its host galaxy. In analogy to the case stud- 
ied by Witt & Gordon ( 1996|, this could provide a good descrip- 
tion of the structure and evolution of the absorbing medium and 
help constrain the opacity and filling factors of the dust geometry 
and clumps from the distribution of Ay™ values in star-forming 
galaxies out to very high redshift. 

The hosts of the dustiest afterglows provide a different pic- 
ture of GRB host galaxies compared to the hosts of optically- 
selected bursts. Although both samples are overlapping in their 
properties, the galaxies in this work have typical luminosities of 
around L* and stellar masses of ~ IO^^Mq, more luminous 
and massive when compared to the hitherto discovered faint and 
blue hosts. Although the number counts are still low, this work 
indicates that a selection based on a large Ay™ picks up pref- 
erentially the more massive and chemically-evolved GRB hosts, 
which is in qualitative agreement with searches for dark GRB 
hosts fP erley et al.|2010a|[20TTa} [Rossi et al.|201 Ih) . 

This suggests that the properties of complete GRB host sam- 
ples are diverse, and complex selection biases are still present: 
not only are the very faintest GRB hosts missing due to inherent 
sensitivity limits, but also some of the brightest, most luminous, 
and chemically evolved ones, because they have not been lo- 
calized accurately enough. Fairly large and massive, dusty, and 
metal-rich galaxies are able to host GRBs, and the trend of low- 
metallicity GRB hosts is not as significant as claimed in previous 
studies, and possibly a selection effect of the young galaxy pop- 
ulation dominating the global SFR at low-z (e.g., Berger et al. 
|2007[[Mannucci et al.|2010b ). This has substantial implications 
for the feasibiUty of tracing the star-formation history with GRB 
hosts and also for the progenitor channels of GRB production as 
a result of their metallicity dependence. In the former case, this 
work indeed indicates that the deficiency of high-mass GRB host 
galaxies in previous studies was at least partially due to a selec- 
tion bias. The latter case, however, depends quite strongly on the 
assumption that the host-inferred metallicities are representative 
of the composition of the progenitor star, while different sight- 
lines through a GRB host can show a dispersion in metallicity of 
around a factor 100 ( [Pontzen et al.|2010| l. 

Intriguingly, all GRBs with Ay™ > 4 mag have metals-to- 
dust ratios significantly below what is typically measured for 
GRB afterglows, and more in line with measurements from the 
Local Group. In addition, there is a strong anti-correlation be- 
tween the metals-to-dust ratio and the visual extinction along 
the GRB sight-line. This effect seems independent on the spe- 
cific host properties and can be interpreted as evidence of two 
physically independent absorbers: dust-free, ionized metals in 



the circumburst environment (and probed by soft X-ray absorp- 
tion), and in contrast a dusty absorber further out (and probed by 
reddening measurements in the UV/optical/NlR). 

A dust column independent of the immediate circumburst 
environment is further supported by the relation between af- 
terglow Ay™ and host properties, in particular the on-average 
higher stellar mass and redder colors. Coupled with the blue 
and very young hosts of, e.g., GRBs 080605 or 100621A it in- 
dicates a complex dust geometry with different constituents in 
the diffuse ISM and in localized patches, which are plausibly 
associated with the cold-neutral absorber detected in rest-frame 
UV/optical spectra. 

Further advances can now be made by getting direct observa- 
tional access to more dusty sight-lines including Ay™ > 10 mag 
events at an increasing redshift interval. Similar observations for 
a large enough sample would investigate the dependence of the 
global dust enrichment with cosmic evolution, and constrain the 
fraction of dust-enshrouded star-formation out to very high red- 
shifts. A sophisticated observational strategy coupled with state- 
of-the-art instrumentation makes such a challenging study fea- 
sible: A rapid response of the order of several minutes by a 
NIR imager at an 8m-class telescope would have enabled the 
detection of the afterglows of GRB 100621A (z ~ 0.5) up to 
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30 mag, GRB 081109 (z ~ 1) up to A'^™ ^ 20 mag, 
and GRB 070802 (z ~ 2.5) up to A'^™ ^ 10 mag, which subse- 
quently could have been followed-up using NIR spectroscopy. 
Once an accurate localization as well as detailed information 
about the GRB sight-line is available, their hosts should be read- 
ily accessible for multi-wavelength surveys via large ground- 
and space-based facilities, yielding information about otherwise 
fully extinguished environments and unprecedented insights into 
the conditions of star-forming galaxies throughout the Universe. 
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